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Abstract. We use the recently derived form factor expansions of the diagonal 
two-point correlation function of the square Ising model to study the susceptibility 
for a magnetic field applied only to one diagonal of the lattice, for the isotropic 
Ising model. We exactly evaluate the one and two particle contributions x S an d 

(2) rr 

X d of the corresponding susceptibility, and obtain linear differential equations for 
the three and four particle contributions, as well as the five particle contribution 
xj, 5 ' ('-), but only modulo a given prime. We use these exact linear differential 
equations to show that, not only the russian-doll structure, but also the direct 
sum structure on the linear differential operators for the n-particle contributions 
(n) 

\ d are quite directly inherited from the direct sum structure on the form 

factors We show that the n th particle contributions x d nave their 

singularities at roots of unity. These singularities become dense on the unit circle 
| sinh2_B„/fcTsinh2i?, l /fcT| = 1 as n -» oo. 



PACS: 02.30.Hq, 02.30.Gp, 02.30.-f, 02.40.Re, 05.50.+q, 05.10.-a, 04.20.Jb 

AMS Classification scheme numbers: 33E17, 33E05, 33Cxx, 33Dxx, 14Exx, 
14Hxx, 34M55, 47E05, 34Lxx, 34Mxx, 14Kxx 

Keywords: susceptibility of the isotropic square Ising model, two-point 
correlation functions of the Ising model, singularities of the square Ising model, natural 
boundary, Fuchsian linear differential equations, complete elliptic integrals. 

1. Introduction 

The magnetic susceptibility \ of the two dimensional Ising model is expressed in terms 
of the two point correlation function C{M, N) — (ao,oCM,./v) as 



oo 




(1) 



M=-oo N=-oo 
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where M(0) is the spontaneous magnetization (which vanishes for T > T c ). This 
susceptibility has been studied for over three decades by use of the form factor 
representations of the correlation function [T]-[5] 

C(M,N)± = {1-t) 1 / 4 -J2 cin) ( M ,N) (2) 

n 

where the subscript +(-) denotes T > T c (T < T c ). For T > T c the variable t is t = 
(smh2E v /kTsmh2E h /kT) 2 and for T < T c is t = {sinh2E v / kT smh2E h / kT)~ 2 
where E v and Eh are the vertical and horizontal interaction constants. The sum is 
over odd (even) values of the integer n for T > T c (T < T c ) and the C (n) (M, N) are 
explicit n fold integrals. Using this form factor decomposition ([2]) in the expression 
for the susceptibility (TTJ we find 

kT- X ± = (l-i) 1/4 'Ex W (3) 

n 

where the x can be expressed in terms of double sum of the form factors: 

00 00 

X in) = E E C^{M,N) (4) 

N=-oo M=— 00 

The study of these x was initiated in [T] with the explicit evaluation of the 
integrals for x an d X ■ However, a real understanding of the analytic structure 
of x began only in 1999 and 2000, with the demonstration by Nickel [HE], for the 
isotropic Ising model (E v = Eh = E), that x^ has a set of singularities, lying on 
the unit circle | s'mh2E /kT\ = 1, which become dense in the limit n — * 00. If these 
singularities do not cancel in the full sum ^ then the susceptibility will have a natural 
boundary at | sinh2£ , /fcT| = 1. 

The existence of a natural boundary in the susceptibility is a profound effect 
not envisioned in the traditional scaling, and renormalization, description of critical 
phenomena. In order to obtain further insight into existence of such natural boundary, 
several of the present authors have made a detailed study of and x^ in [6]- [9]. 

The integrals for x , as explicitly written out in [H[S], are quite complicated, 
and, for that reason, it is difficult to extend the analysis of [6]-[9] to x^ with n > 5. 
A direct attack [10] on x^ indicates that more than 6000 terms and, for x^\ probably 
more than 20000 terms in the power series expansion in t are needed in order to find the 
linear differential equations which they satisfy. Therefore it would be most useful to 
study "model integrals" that are simpler to analyze, and which incorporate significant 
features of x^ ■ Several such "model integrals" have been previously studied [11] and 
revealed a rich structure of singularities beyond those found by Nickel [H [5] . 

In this paper we introduce what is probably the most physical simplification 
of the Ising susceptibility which retains the property of having singularities on the 
unit circle |i| = 1. This model is obtained by considering the isotropic Ising model 
with a magnetic field which acts only on one diagonal of the lattice. The magnetic 
susceptibility for this diagonal field will be, then, given by the diagonal analogue of 
Q: 

oc 

kT 'Xd = E ( C ( N ' N ) -M(0) 2 ). (5) 
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This expression^! should be much easier to study than the full susceptibility. This 
stems from the form factor decomposition of the diagonal two-point correlations 
C (N, N) , that has been, recently, presented [T2] and proven [T3] , and which is much 
simpler than the decomposition obtained directly from [T]. In particular for T < T c 



C(N, JV)_ = (1 - t) 1 / 4 • (l + f (2n) ( N > *)) ( 6 ) 



with 



t n(n+\N\) l 



1 2n 



X 

3 



n 



(n!) 2 7T 

X 2 j-l(l-X 2 j)(l-te 2 j) 



5fe • x 



k 



k=l 

1/2 



i V^jCl — rc 2j -i)(l - tx 2] -i), 

x n n ( j - tx 2j-ix 2k )- 2 (7) 

l<j<n l<k<n 

X II ( X 2j-1 - X 2 k-l) 2 (X2j - X 2 k) 2 
l<j<k<n 



and for T > T r 



C(N,N) + = (l-t) 1/4 -X;/ (aB+1) W*) ( 8 ) 

n=0 

with 

£Ti(n+l)+|iV|(n+l/2) ,-1 2n+l 



f (2n+l) {Nt) = 1 / ... / TT ^.JAI 

^ v ; 7r 2 ™ +1 n!(n + 



hl/2) /-l 2n+l 

if ■/„■■■ / n 



X 



JJ ((1 - »2j)(l ~tx 2 j) ■ X 2 j 



1/2 



X 



i=i 

TT / x i 2 



x n n c a -^aj-i^fe) -2 (9) 

l<j<n+l l<fe<n 

X J| (^2j-l - ^2fc-l) 2 • JJ (X 2 j ~ X 2k ) 2 

l<j<k<n+l l<j<k<n 

From now on, the integer N should be understood as \N\ when evaluated. 
Thus, if we use ^ and ([5]) in ©, and evaluate the sum on N as a geometric 
series, we obtain for T <T C 

oo 

kT- X< l-{t) = (l-i) 1/4 -Exl 2n) W (10) 

n=l 

with 

x y d '(*) 



(nl) 2 ^ 2n Jo Jo M 1 -t^i---x 2 „ 



If Such partial sums have already been introduced and the asymptotic expansion of their short term 
part considered in Orrick et al. I14| . 
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X 

3 



*2i-l(l - X2j){l ~ tX 2 j) N ' 



X 2 j(l - X 2 j-l)(l - tx 23 -l), 

n n ^ ~~ tx 2j-\ x2kY 2 

l<j<n l<k<n 

Yl (X2j-1 - X 2 k-l) 2 (X2j - %2k) 2 (11) 
l<j<k<n 



and for T > T c 

kT- Xd+ (t) = (l-t) 1/4 -£^ 2 " +1) W (12) 

n=0 



with 



1 2n+l 



j.n(n+l)) /"l 

W = ^n\(n+iy.-J -L 11^* 



1 +t"+ 1 / 2 x 1 --- a ; 2Tt+1 " / X i/2 
l-^+V2 Xl ... X2 „ +1 ' II ((! - ^0(1 - t^) • ^ 



j=l 



-l 

< 1 ~ 

i=i 

x n n ( x -^■-^r 5 

l<j<n+l l<fe<n 

X JJ (X2 3 -1 ~ X 2 k-l) 2 Yl ( X 2j- X 2k) 2 - 

l<j<k<n+l l<j<k<n 

The expressions and (|13[) are, indeed, much simpler than the corresponding 
expressions for x given in 

In this paper we analyze the contributions Xd M ^° the diagonal susceptibility. 
In section ([2]) we analytically evaluate the integrals for n = 1 and 2. In section ([3]) we 
present, and analyze, the ordinary linear differential equations satisfied by x d W an d 

(*)) as we ^ as Xd\t) modulo a given prime. These linear differential equations 
have a direct sum decomposition which we relate in section ([5]) to the direct sum 
decomposition found in [T3] for the form factors /( n )(7V, i). In section (fj| we extract 
the singularities in the integral representations (fTTj) and (fTBf for -)Q it) ■ We find that 
all the singularities of (t) are at the roots of unity t n = 1, while, for x^ 2 ™ +1 ^(£), 
they are at i™+ 1 / 2 = 1. These singularities for t ^ 1 are the counterparts, for Xd^i^), 
of the singularities on the unit circle found by Nickel [H[S] for the full susceptibility X- 
However, unlike X (t) which was shown in [5]-[llj to have many other singularities 
which lie outside the unit circle | smh2E/kT\ = 1, the diagonal Xrf^W nas no further 
singularities other than those which satisfy t n = 1 or i n+1 / 2 = 1. We conclude 
in section © with a discussion of the several different types of singularities in Xd(t) 
which follow from the singularities in Xj^\t)- 
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2. Evaluation of x ( p{t) and xf\t) 

The contribution of Xd\t) is explicitly given from (fT3"|) as 

X ( a\t)= ^ j\ x .\±^±.[{l- x ){l-tx)x]-^ (14) 

which, setting x = sin 2 is written as 

m, N 2 f^ 2 1 +i 1 / 2 sin 2 6» , , N , /9 

Recalling the definitions of the elliptic integrals of the first kind 

= / n i ■ 2 flU /2 ( 16 ) 

Jo (1 - i sin Oy' 2 

and of the third kind 

ni( "'* )= X (l+,sin 2 0)(l-;sin 2 e)V 2 
we may write (fl"5|) as 

xi 1} (t) = ^•(2n 1 (-t 1 /a,t) -A-(i)). (18) 

Thus, if we use the identity [T5] 

- • (lix{v, t) +ni(t/i/, t)) = (19) 



7T 



with j/ = — t 1 / 2 we find: 



x { d 1] (t)= T-^r/a- ( 20 ) 

It is instructive to derive (|20|) without recourse to identities on elliptic integrals. 
We first rewrite (|14|) as a contour integral 

*£ )(t) = 2^ / f ■ T=f ■ [(1 _ * 1/2z)(1 _ <1/2 ^ 1 )] _1/2 ( 21 ) 

on the contour |z| < 1 (which is, in fact, the form in which the form factors f( n \N, t) 
were originally derived |13j). The integrand in (|2ip is antisymmetric if we send 
z — > 1/z. Therefore 

X^W = ~xi 1} W -(residue at z = 1) (22) 
The residue at z = 1 is easily evaluated as 

2 

(residue at z = 1) = (23) 



and thus using (f23|) in (f22|) we again obtain the result (f20|) . 

To compute Xd (*) we a ls° use the contour integral method, and rewrite (fTTj) 
with n = 1 as 

xi 2) W = j^fdztjdz* [24 1 

1/2 



1 + Z±Z 2 
X (1 - Z!Z 2 ) 2 



(1 -t l / 2 z 2 )(l -t ll2 z^) 



(1 -tV2 Zl )(i -tVa^j" 1 ) 
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on the contour \z^\ < 1 for k = 1,2. As was the case for x*p{t), we note that the 
integrand of (I24[) is antisymmetric if we send z\ — > l/z%, z 2 —* 1/^2- Therefore 

xTi 1 ) = ~ (2^)2 f ^i( residue at z 2 = Zi 1 ) (25) 

where: 

(residue at z 2 = zf 1 ) = zf 3 [(1 - < 1/2 zi)(l - i 1 / 2 ^ 1 ]" 1 / 2 (26) 

X ^ {(1 + ZlZ2)[(1 -* 1/2 ^)(1 -* 1/2 ^ 1 ] 1/2 }l„- 

When (|26|) is evaluated, and substituted into |25|) . the resulting integral over z\ has 
only poles. Keeping only those terms which give nonvanishing contributions, we find 

*?(*) = ±t dzi —-—J—-— = -J—. (27) 



8iriJ 1 (1 -tV*zi)(zi -t 1 / 2 ) 4(l-t)* 

It should be noted that neither x rf (£), nor Xd contain logarithms even though 
there are logarithms in both / W (JV, i) and /( 2 ) (AT, i) for all N. This is to be contrasted 
with the corresponding results for the full susceptibility where it was seen in [T] that 
jC n '{t) has no logarithms for n = 1, but does have a term in hit for n = 2. 

3. Linear differential equations for PQ (i), Xd\t) an( i Xd W 

We now turn to Xj(t) for n > 3. When written in contour integral form 
the integrands still have the property of being antisymmetric when z^ — > 1/zk- 
Unfortunately this property is not sufficient to reduce the computation to integrals 
that can all be evaluated by residues. Therefore we do not have an explicit evaluation 
in terms of elementary functions, and we continue our study by using formal computer 
computations to determine the linear differential equations satisfied by Xd W > Xd (*) 
and Xrf 5 '' (*) ' as was done for x^ an d X^ f° r the full susceptibility in [6]-[9] and for 
the (N, t) in [12] . We present the results of these computer calculations for X4 (*) > 
X { d ] (*) and Xd } (*) separately. 

3.1. Linear differential equation for x d {t) 

For Xd^ we chose x — t 1 ! 2 — smh.2E v /kT smh2E} l /kT as our independent variable. 
We find that the linear differential operator for x d ( x ) i s °f order six, and has the 
direct sum decomposition 

4 3 > = 4 3 >e4 3) e4 3) (28) 

with 

L\ 3) = Dx + — , (29) 

35—1 

L f = Dx * + 2 ( 1 + 2a; ) . Dx + , (30) 

2 (l + a;)(x-l) (1 + at) (a;-l)a;' V ; 
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4 3) = Dx 3 (31) 



^3 

3 (8 x 6 + 36 x 5 + 63 a; 4 + 62 x 3 + 21 x 2 - 6 x - 4) 
+ 2 (a; + 2) (1 + 2 x) (1 + x) (x - 1) (1 + x + x 2 ) x 

+ -! 

(x + 2) (1 + 2 a;) (1 + x) 2 (x - I) 2 (1 + x + x 2 ) x 2 

no 

[x + 2) (1 + 2 a;) (a: - l) 3 (1 + x + x 2 ) (1 + x) 2 x 2 



■ Dx 2 
Dx 



with: 



no = 2 a; 8 + 8 x 7 - 7 x 6 - 13 a; 5 - 58 x i - 88 a; 3 - 52 a; 2 - 13 a; + 5, 
ni = 14 a; 8 + 71 a; 7 + 146 a; 6 + 170 a; 5 + 38 a; 4 
- 112a- 3 - 94x 2 - 19x + 2. 
The solution of is in fact (up to a constant) Xd\t) as given in (|20|) . This 
means that x d (t) is actually solution of the Fuchsian linear differential operator >Q 

corresponding to \ d (t) . We recover then the "russian-doll structure" noticed in [8l [9] 
for the third contribution x^ to the full susceptibility. 

(3) 

The linear differential operator of order two, L\ has, at x = 0, 1, — 1 and x = oo, 
respectively the exponents (0, 1), (—2, 0), (0, 0), and (1, 2). Some manipulations on the 
formal solutions, at x = 0, give the result that the solution analytic at x — reads 



sol{Lt>) = — • K(x 2 ) + • £(^) (32) 

with: 

K(y)= F (1/2, 1/2; 1; j/), = F (1/2, -1/2; 1; j/) (33) 

(where we have abused conventional notation by omitting the factor of tt/2 in the 
canonical definition^ of the complete elliptic integrals K and E) . 

f 3) 

We have not found the explicit solution of the linear differential operator L y 3 ' 
which has the following regular singular points and exponents: 

l + .T + a; 2 =0, p = 0,1,7/2 -> a; 7/2 , 

a: = p = 0,0,0 — > log 2 terms, 

x= 1 p = -2,-l,l-> a;" 2 , a;" 1 , (34) 

x — —1 (0 = 0,0,0 — > Zog 2 terms, 

x = oo (0=1,1,1 — » Zog 2 terms. 

The singularities at x = 2, —1/2 are apparent. 

5. Linear differential equation for \^ (t) 

For (t) we use t as the independent variable. The linear differential operator for 
X d f (t) is of order eight, and has the direct sum decomposition 

4 4) = 4 4) ©4 4) ©Lf (35) 

I In maple's notations, for t = k 2 (k is the modulus), K(t) = K(k 2 ) in (133 I t reads : 
/ij/pergeom([l/2, 1/2], [l],t) = 2/-7T ■ EUipticK(k), but reads 2/7r ■ EllipticK[k 2 ] in Mathematica. 
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with 



^=nt +¥ - Wt , (36) 



(5t 2 + 6t-l) n+2 (3t 3 + 6t 2 -2t-l) 
_ 

3 



d 4} = Dt A + ^- • Dt z + ±- — ^- ■ Dt 

3 (1+t) (*-l)* (l + t)t 2 (t-l) 2 



2(l + t) (i — 1) £ 2 ' 

m d (7i 4 -68i 3 - 114t 2 + 52i-5) , 

4 + 1) (t 2 - 10 i + 1) (t - 1) i 

C5 1 5 - 55i 4 - 169 i 3 + 149 1 2 - 28t + 



(37) 



{t + i){t 2 -iot + i)t 2 (t-r 



2 



Z?i 2 



(i 4 - 13t 3 - 129i 2 + 49i - 4) , 

I 2 -i 4r • 38) 

(i + 1) (i 2 - 10 £ + 1) £ 2 (t - l) 2 



(£ 2 - lOt + 1) (£ ~ l) 2 £ 3 ' 



The solution of is (up to a constant) the function x)p(t) given in (|2T|) . Here 
also, (f ) is solution of the Fuchsian linear differential operator corresponding 

to X^cP{t)-> an d we recover, again, the russian-doll structure noticed in [51 for the 
fourth contribution x^ to the full susceptibility. 

The solution of L3 is found to be (with notations (133)) ): 

S0 Z(4 4) ) = -K(t) 2 + j±±L . E {tf + • K(t) E{t). (39) 

We have not found the explicit solution for the operator which has the 

following regular singular points and exponents 

t= 0, p = 0,0, 0,1 -> log 3 terms, 

t= 1, p=— 2,— 1,0,1 — > i~ 2 , t _1 , fog term, 

* = -1, (0 = 0,1,2,7 -> t 7 log term, (40) 

t = oo, p = 0,0, 0,1 — > iog 3 terms. 

The singularities at the roots of i 2 — 10< + 1 = are apparent. 

3.3. Linear differential equation for xf' \t) 

For Xd^ we cnose : again, x = t 1 ! 2 — smh2E v /kT smh2Eh/kT as our independent 

variable. The first terms of the series expansion of xf^ i x ) rea d : 

-("51/ n 3 i, 39 14 5085 lfi . 

vy ; (a;) = -a; 12 H x 14 H -a; 16 (41) 

Ad v ; 262144 1048576 67108864 v ' 

9 ^ 33405 _ il8 , 315 ^ , 



67108864 268435456 536870912 

Now, the calculations, in order to get the linear differential operator for x^\ x )i 
become really huge. For that reason, we introduce a "modular" strategy which 
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amounts to generating large series modulo a given prime, and then deduce, from 
a Pade-Hermite procedure, the linear differential operator for xf^ ( x ) modulo that 
prime. We have generated 3000 coefficients for the series expansion of X d ( x ) modulo 
a given prime (here 32003), and actually found a linear differential equation modulo 
that prime of order 25, 26, • • •, using 2200 terms in the series expansion. One can 
also obtain linear differential equations of smaller order for but where more 

terms (2500, 2600, 2800, • • •), are needed. The polynomial corresponding to apparent 
singularities in front of the highest derivative is now very large. The smallest order one 
can reach is 19, and the linear differential equation we have "guessed" has required 
more than 3000 terms in the series expansion. Note that a linear differential equation 
of the minimal order is not, necessarily, the simplest one, as far as the number of terms 
in the series expansion needed to guess it, is concerned. We have already encountered 
such a situation in [9]. 

Recalling C & , the order six linear differential operator corresponding to \ d \ x \ 
one finally finds that the linear differential operators we have obtainedj for 

Xd (x), have the following factorization: 

4 5) = 4fV4 3) (42) 

implying that Xd ( x ) ^ s ac t uau y a solution of C„ \ the linear differential operator 
for x d 5) ( x )- Thc 

"russian-doll" structure shows up again. For the linear differential 
operator Cn of smallest order (n = 19), 4-6 nas a polynomial of apparent 

(5) 

singularities of degree 331. For larger orders for C n _ 6 we get smaller apparent 
polynomials. In |Appcndix A[ we sketch an order twenty linear differential operator 
4 5 -6 modulo the prime 32003, which has no apparent singularities, and only thc 
"true" singularities of the linear differential operator Cn" 1 for Xd ( x )'- 

(x + 1) • (x - 1) • (x 2 + x + 1) • (x* + x 3 + x 2 + x + 1) -x = 0. 

(5) 

We have not yet been able to get the direct sum decomposition of £„ ' in this approach 
modulo prime. 

4. Singularities of x|^(i) 

In [12] we found that the form factors f( n \N, t) have singularities of the form ln n (l— t) 
at t — > 1. These singularities come from the factors [(1 — Xj)(l — txj)] 1 ^ 2 in the 
integrands of (J7J) and These factors are also present in the integrands (TlTj) and 
P3[) for X d n \t)t an dj thus, we expect that, for general values of n, there will be powers 
of ln(l — t) present in x d n ^ e see ' fr om the previous section, that these logarithmic 

('X) 

singularities are first seen in the linear differential equation for x d (£)• 

However, there are additional singularities in Xd^ (*) coming from the vanishing 
of the denominator 

1 - t n xix 2 ■■■ x 2n , (43) 
and, in Xd +1 C0i coming from the vanishing of the denominator 

1 - t n+1 / 2 x lX2 ■■■x 2n+1 (44) 

| The method consists in searching, modulo a given prime, the linear differential operator for 

4 3) (xi 5) ). 
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which are not present in f^ n \t). For X? the vanishing of ([4"5)) occurs for t n = 1 
at the endpoints X\ = x 2 = ■ ■ ■ = x 2n = 1, and, for ^ 2 " +1 \ the vanishing of (|44p 
occurs for t'^ 1 / 2 = 1 at the endpoints x% = x 2 — ■ ■ ■ = x 2n +i = 1. 

When t — » 1 this additional singularity is the simple pole (1 — for both 
Xd (t) an d Xd 2 " +1 ^(*)- When £ approaches the roots of unity t" = 1 (with t ^ 1) 

t JjT , = e 27ri V« w ith Z = 1, 2, Ti-l, (45) 
then x^ 2 ™' CO nas a singularity of the form 

«2„- (*-ti,„) 2 " 2_1 -ln(t-^„). (46) 

Similarly, when £ approaches the roots of unity for x^ 2 ™^ (with ( ^ 1), namely 
t[] +1/2 = 1, then x (2n+1 \t) has a singularity of the form n 2n+1 ■ (t 1 / 2 -tJ/ 2 )(»+i) a -V2. 

These singularities, on the unit circle in the complex t plane (for t ^ 1), are the 
counterparts for the diagonal susceptibility of the singularities found by Nickel [H [5] 
on the unit circle \s\ — 1, with s = sinh2i?/fcT, for the x^ of the isotropic Ising 
model. 

There are no other values of t in the complex plane for which the functions xk W 
are singular. This is in distinct contrast with the x ( s ) °f the isotropic Ising model 
which have singularities at many other places on the complex s plane [B]-[S],[TT]. In 
this section we derive the behavior of X^(i) & ^ ^ e ro °t °f unr ty points t n = 1 and 
f n+i/2 _ L Appendix B gives the values of the amplitudes at the singular points lying 
on the unit circle \t\ = 1, for jffl and Xd > usm S the matrix connection method [§]• 

4-.1. The singularity in x^ at t = \ 

To extract the dominant singularity in x^ 2 ™^ at i = 1 we set t = 1 — e and 
Xk — 1 — e ■ yk in and set e = wherever possible. Thus we obtain the 

result that as t — > 1 (e — > 0) 

XT" ~ (47) 



with: 



J f = ^2^2n / ^1 ' ' ' / ^ ■ I] 



V2j(i + y2j) \ 1/2 



— n n (i+w«-i+itt*) 

n + yi H hj/2n . , 

\<]<n l<fc<ri 



X 



-2 



X II {y2j-i~y2 k -i) 2 {v2 ~v2k) 2 . (48) 

l<j<k<n 

~(2n+l) 



Similarly, we find the singularity in Xd ™ at t = 1 reads 



with: 



xT +1) ~ Y~ t • 4 2n+1) (49) 



2 /"OO roo 

n!(n + l) W «Wo ^ ■■•./„ 
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1 n n+1 

n +yi + ... +y2n m^y*)y^ 2 Ife- 1 (i+^- 1 )r 1/2 

x n n ( i+ y2,-i +v2kT 2 

l<j'<n+l l<fc<n 

x n (y^- 1 - y^if ■ n (y 2j -2/ 2 fc) 2 - (50) 

l<j<fc<n+l l<j<fc<n 

^.,2. Singularities for x d 2n \t) at t n — 1 /or t ^ 1 

When n = 2 the root of unity singularities <™ = 1, which is not f = 1, occurs at t = —1 
where the analysis of the differential equation given in section ([3]) shows that there is 
a singularity of the form (1 + t) 7 ■ ln(l + t), (see Appendix B). 

To demonstrate, for general n, that the singularity in X d 2n \t), at the points 
t n = 1, is given by (|46|) . we set £ = ■ (1 — e) and Xk = 1 — e ■ 2/fe. Furthermore, 
because the singularity only occurs in a high derivative of X^ 2 ™" 1 ^), we consider the 
m th derivative of X^ 2 ™^), and will eventually see that m should be chosen to be 
2 n 2 — 1. Then using i = t; iTl • (1 — e) and = 1 — e • in the m th derivative of 
(jlip . and setting t = i; iTl and Xk = 1 wherever possible, we obtain a model integral 
whose leading singularity, at t = ti m , will be the same as the leading singularity in 

the m th derivative of x d 2n \t): 

,2n 2 -l-mo m | foo roo » 



(n!)2 tt2« (1 - ti, n ) 

X (^i)V2 . 



J'=l 



V2j (n + yi+ij2 + --- + y2n) m+1 

X II (^2j-l - V2k-l) 2 (y2j ~ Z/2/c) 2 - 
l<j<fe<n 

This expression is formally independent of e when m = 2n 2 — 1. When m has this 
value the integral diverges logarithmically when all y^ become large, such that all 
ratios i/j/yk are of order one. We thus conclude that the singularity in X^ 2 "' ) (£), at 
t = ti !n , is given by (gBJ. 



4.3. Singularities for Xd (*) at t n+1 ' 2 = 1 for t j= 1 

For the function Xd (*) we use d in section ([3]) the independent variable x — i 1 / 2 , 
and found, at the points Xq + xq + 1 = 0, that there is a singularity of the form 
(x — a;o) 7 ^ 2 , (see Appendix B). An analysis, completely analogous to the analysis 
given above for x^ 2 "^); demonstrates that, for all n, the singularity in Xd™ +1 ^(*)> at 
xl n+1 = t l+1/2 = 1 (with t± 1) is given by K2 „+i • (x - x ) {n+1)2 " 1/2 . 



5. The direct sum structure 

Perhaps the most striking feature of the 

Xd 4 \t) exhibited in section (PJ is their 
decomposition has previously been seen 



linear differential operators for x d (t) an d 
decomposition into a direct sum. Such a 
for the full susceptibility where the linear 
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differential operator for x^ is the direct sum [B] of the differential operator for x' 1 ) 
and a second linear differential operator. Similarly, the linear operator for x^ is the 
direct sum [H] of the linear differential operator for x' 2 ' and a second linear differential 
operator. In the papers [B]-[5] the question is posed of how general is the phenomenon 
of the direct sum decomposition. 

In our previous paper [12] on the form factors (A, t) of the diagonal Ising 
correlations we found that the linear differential operators of all form factors for n < 9 
have a direct sum decomposition, and that this decomposition is surely valid for all 
values of n. 

An inspection of the direct sum decompositions reveals a great deal of structure 
which is relevant to yffi (t). As an example consider the odd form factors /( 2n + 1 ) (TV, t) 
relevant for T > T c . From the direct sum decomposition of the linear differential 
operators [T2] for /( 2,l+1 )(A, t), we find the following form: 

fW(N,t) = + f (1) (N,t) +g (3 \N,t), (52) 



. 2 6 . 

f {5) (N,t)= ± ■ (15 A 2 + 40 A + 9) ■ /« (A, t) 

+ i(A + l). 5 ( 3 '(A,i) + g^(N,t) (53) 

where t N / 2 -gW(N,t) is a homogeneous polynomial of degree n in the complete elliptic 
integrals K and E, with coefficients which are polynomials in t. If we then sum the 
decomposition flS^) of fW(N,t), we find that 

X ( i\t) = +x { ll(t) +m (54) 



wherd 



*S(*) - \ E = \&\ (55) 

JV=— OO 

°° AT 

x%= E (56) 

7V=-oo 

where \ 1& * ne solution of the linear differential operator (which is regular at 
x = t 1 / 2 = 0) and where 



~(3) 



9 {3) (N,t) = h-x^(t) +sol(L ( i ) ) (57) 

JV=-oo 

witri bi ^ -1/6. Thus we see that the direct sum decomposition of the 
linear differential operators for x d (t) follows immediately from the direct sum 
decomposition of the differential operators for /( 3 )(A, t). 

For /( s )(JV,t), the sum (|53]) can be written in terms of f^{N,t) as 

f (5 \N,t) = ~ (15 A 2 +l)-f^(N,t) 

+ l(A + l)./( 3 )(A,t) +g^(N,t) 

% The integer TV should be understood as \N\ in all the summations below. 

(I In fact, and since the full factorization of is known, the projection of jq [i) in xS (*) can b c 
computed and one finds fei = 1/6. 
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Summing on N, one obtains: 



(5) 
Xd 



N N N 

lj2N-f (3) (N,t) +J29 (5) (N 7 t). (58) 



N N 

The first term, at the right-hand-side of (|58|) . satisfies the linear differential equation 
of order-one corresponding to Xd \ (i- e - m (EHJ) ) - The second term satisfies an 

(3) 

order-one linear differential equation. The third term is just , up to a constant, a 



result that we knew from the modulo prime method of section (|3.3|) . We have found 
the differential operator corresponding to the fourth term which is of order eight and 
has, in its direct sum, the linear differential operator given in (|28|) . As was 

the case for the sum in (|57|1 . one might imagine that the linear differential operator 
corresponding to the last term in ([55]) will contain some of the differential operators 
related to the previous terms. 

We consider now the x^ 2 "^)- In [12] we found that there is a direct sum 
decomposition for the linear differential operator for /( 2 ™) (N,t), just as there was 
one for /( 2n+1 )(iV, t). From these direct sum decompositions, we find: 

f^(N,t) = ^ +g (2) (N,t), (59) 



f*>(N,t) = ^-^ + U + V 9 (N,t) +9 (N,t) 

= ^^ + ^^)-f^N, t) + 9 ^N, t) (60) 

Similarlji, the direct sum decomposition of the linear differential operators for 
f {2n \N,t) will lead to the direct sum decomposition of the differential operators for 
X d 2n \t) for all ri. 

We conclude that there will be a direct sum decomposition for the linear 
differential operators for Xj^ft), for all n, and that this direct sum decomposition 
is inherited from the direct sum decomposition for the linear differential operators of 
the form factors of the Ising model. 

Remark: It is worth noting that we have also performed a large set of calculations 
(that will not be detailed here) on more "artificial" toy susceptibilities like : 

oo oo 

Xtoy = EE N*.f&(N,t). (61) 

3=1 N=l 

Similarly, for the corresponding j-particle contributions Xtoy we f° un d for the 
first values of j (j = 1, 2, 3, 4), the corresponding Fuchsian linear differential 
operators. The singularities of these Fuchsian linear differential equations and of 
the corresponding j-fold integrals, are totally and utterly similar to the one's of the 
diagonal susceptibility analyzed in this paper. Again we have equations totally similar 
to (|54p. (fS"5")) . (fS"6"]) . This confirms, very clearly, that the direct sum decomposition for 
the j-particle contributions Xtxyy ^ s straightforwardly inherited from the direct sum 
decomposition for the form factors of the model. 

X Note that the first terms in Il59|l or l|60|l do not lead to divergent sums. They are balanced by the 
sums on the last terms, i.e., g^{N,t) and gM(jV, t). 
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5.1. Resummations 

There is one further feature of these decompositions which must be mentioned. 
Namely, the sums of the form 

oo 



^ N p ■ f^ n+L >(N,t) (62) 



^ N p ■ /( 2n+1 

N=-oo 

(3) 



will diverge at t = 1 as (1 — t 1 ^ 2 ) p 1 . Thus, for example, the solution of L 2 , given in 
(f32jl . diverges at t 1/2 = 1 as (1— t 1 ! 2 )^ 2 . We see, from ([34]) . that the leading singularity 
in the solution to Lg will also diverge, at t 1 ! 2 = 1, as (1 — £ 1//2 )~ 2 . However, the full 

solution for (i) must diverge, at t — > 1, only as (1 — t 1 / 2 ) -1 , and, therefore, there 
must be cancellations between terms in the direct sum. This phenomenon will happen 
for all X^ n+1 \t) where the cancellations become more extensive as n increases. For 
heuristic reasons let us consider Xg versus f( l \N). 

The diagonal susceptibility of order one is defined as 

oo oo 

xf = E f W (N,t) = fW(Q,t) + 2 V fM(N,t) (63) 



/ . J V ' / J V > 7 1 / . 

^=-00 JV=1 



with 



f W (iVl t) = ■ t N ' 2 ■ 2 F, (1/2, 1/2 + N; 1 + N; t) . (64) 

Writing in the sum over N, the hypergeometric function as a series, one obtains: 

k=ON=l K )k ■ ■ 

with (x)k denoting the Pochhammer symbol, and K{t) is as defined in (|33[) . One may 
shift N, separate the even from the odd N, and try to sum. One may also generate 
the expansion, find the linear ODE and solve. Remarkably the sum in (f63|) reduces to 
a simple expression in terms of the complete elliptic integral of the first kind: 

2 f; /^(iV,t) = _-^L-J- -K(t). (66) 

N=l 

This remarkable identity (|66p explains why a sum like (|63p . where each term is 
polynomial expression of the complete elliptic integral of the first kind and of the 
second kind, succeeds to reduce to a simple rational expression in t 1 ^ 2 : 

X?'= A'W-j^-^-A-W . _L_. (67) 



6. The singularities in the diagonal susceptibility 

Thus far we have discussed the n th particle form factor contribution of Xatf)- ^ 
remains to use this information to study the diagonal susceptibility Xd(t) itself, and, 
for this, we need to consider several problems which also occur for the computation 
of the full susceptibility x( s )- 

We computed the susceptibilities Xd± by summing the form factor expansions ([6]) 
and ([5]) of the diagonal correlation functions C(N,N) over all integer values of N. 
In doing this we have interchanged the sum over position N with the sum over form 
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factors n. In field theory language we have interchanged the high energy limit with 
the sum over n particle intermediate states. This interchange is universally done in 
both statistical mechanics and in field theory, but should, in principle, be justified. 

Let us assume that this interchange of the sum over N and n can be made. Then 
the behavior of Xd, as T T c , can be studied from the behavior of X^(t) as i — > 1 
if we make the additional assumption that the limit T — > T c can also be interchanged 
with the sum over n. We then may use (|47| in (fT0|) to find, as T — > T c — , that 

oo 

k T - Xd _ ~ (i-r 3/4 -E^ 2n) ( 68 ) 

n=l 

and, similarly, by using (|4T))) in (jT^J) we find, as T — > T c +, that: 

oo 

fcT- Xd+ ~ (l-t)- 3/4 'i;4 2n+1) - (69) 

The sums in (|68| and (|69|) must be shown to converge if these estimates of the critical 
behavior are to be correct. For the full susceptibility, similar convergence has been 
recently demonstrated by Bailey, Borwein and Crandall [To] . 

The single pole divergence, which occurs in Xd W f° r eacn n i 1S the analogue, for 
the diagonal susceptibility, of the double pole divergence (1-s) 2 in the x (n) (s) of the 
full susceptibility. In both cases this divergence occurs from terms in the integrand 
of jQ (t) (or X ^ n \s)) which are not present in the corresponding integrals for the 
form factor representation of the correlation function. These divergences may be said 
to come from long distance effects, and are captured in the scaling theory of the 
correlation functions. 

There are further singularities in yffi (t) which come from the pinch 
singularities [UJ of the square root branchpoints [(1 — x)(l — tx)] 1 ^ 2 in the integrands 
which are also present in the form factors f( n '(N,t). In the form factors f( n '(N,t), 
singularities give divergent terms ln"(l — t) which, term by term, would give the 
dominant contribution (1 - t) 1 / 4 • J2„ hi n (l - *) ^(l - *) to the correlation function. 
However, the overall factor (1 — i) 1 / 4 is absent in the diagonal correlation C(N,N) 
and is thus cancelled by the infinite sum on n. We thus regain the original expansion 
of C(N, N) as an N x N determinant whose singularities, at t = 1, are of the form: 

(1 - t) N * ■ \n N (l - t). (70) 

These singularities in C (N, N) , which come from the summation of the form factors 
over n, may be said to be short distance singularities. 

In x!i (£), f° r n — 1' 2, the results of section @ show that there are no 
logarithmic terms. However, for n = 3, 4, logarithmic terms occur and we presume 
(but have not yet demonstrated) for arbitrary n > 4, that Xd^it) W1 U have logarithmic 
terms. However, just as was the case for fW{N,t), these powers of log's must be 
summed over all n, and, again, just as for f( n \N,t), this sum must cancel the factor 
of (1 — t) 1 / 4 . It will thus give terms of the form (|70p for Xd(t) which are the counterpart 
of the similar terms in the full susceptibility x(s), and have been studied in detail by 
Orrick, Nickel, Guttmann and Perk [14]. 

It remains to discuss the singularities of the diagonal susceptibility on the unit 
circle \t\ = 1. These singularities not only have the property that they become dense on 
the unit circle asn->oo, but since they are roots of unity, they accumulate according 
to a uniform distribution, in contrast with the density of the "nickellian" singularities 
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(see (3.24) in [TJ] and (1) in [IT]), which coincides with the density of zeroes of the 
partition function without magnetic field given in [T7j. Another qualitative difference 
between the diagonal susceptibility, and the full susceptibility, is that these unit circle 
singularities for the diagonal susceptibility are logarithmic for T < T c , and are square 
root type for T > T c , whereas, for the full susceptibility the singularities for T > T c 
are logarithmic, and the singularities, for T < T c , are of square root type. In all 
cases the amplitudes of the singularities depend strongly on n which would seem to 
prevent any possible cancellation between singularities in sets such as ~)^ mn ^{t)-, with 
n fixed and m — - 1, 2, • • •, which have the locations of some of the singularities at 
coinciding positions. Therefore the arguments used by Nickel [l]-[5] to conjecture a 
natural boundary in the full Ising susceptibility will similarly suggest that there is a 
natural boundary in the diagonal Ising susceptibility as well. As is the case with the 
full susceptibility a more rigorous argument would be most desirable. 
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Appendix A. Linear differential operators 

Recalling the factorization relation (|42p of the linear differential operator Cn^ 
corresponding to Xd\ x ) a linear differential operator C^-q an d the linear 

differential operator corresponding to x d ( x ) '■ 



the smallest order for the linear differential operator C n _ & is thirteen, but this yields a 
large set of apparent singularities. For an order twenty, the linear differential operator 

(5) 

C 2Q has no apparent singularities, and requires less terms, in the series expansion, to 

(5) 

be guessed. Let us sketch this order twenty linear differential operator £20 modulo 
the prime 32003: 




r (5) A3) 



(A.1) 



i=20 




(A.2) 



Q20 = (1 + xf {x 2 - l) 3 {x 3 - l) 2 (x 5 -l)x 



(A.3) 



= (1 + x) 9 (x - l) 6 (1 + x + x 2 ) 2 (1 + x + x 2 + x 3 + x 4 ) ■ x 1 
Q m = (1 + aO Q(9 ^ m) ' (x ~ l) a(6 ' m) ■ (1 + x + x 2 ) a{ - 2 ' m ) • x Q ( 10 ' m > • 
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where m — 0, 1, 2, • ■ • , 19 and 

a(N, m) = sup(0, N - 20 + m) 
and where the polynomial^ q m read respectively : 

gig = 29006 ■ (x 2 + 1466 x + 22107) ■ P^ 7) ■ r\ l J ) ■ P§ 3) ■ P\ 



\ The polynomials q m are modulo the prime 32003. 
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q 18 = 23383 ■ (x + 13780) ■ (a; + 13647) • P^ 4) ■ P^ 0) ■ P$ ■ p/ 5 ' 



18 * 18 1 18 1 18 ' 

(59) 



gi7 = 23894 • (x 2 + 13626 x + 7861) • P£ 
9i 6 = 31481 • {x + 31392) • P x (19) • P x (18) • P& 4 > • P^ 
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g 15 = 13122 • (x + 2555) • p£°> • • 

914= 1689 • P(f } • Pi ( 4 10 \ 

913 = 851 • (x + 6919) -P^ 7 \ 

9i 2 = 542 • Pg* ■ Pg* ■ pg 4 > . . p£) . P (3) 5 

9ii = 26141- (x + lJ.P^.pW.pW.P®, 

910 = 31757 ■ x ■ {x + 14054) ■ P§ 5) ■ P$ x > • P x (2) , 

9 9 = 31477 .P 9 (45) .P 9 (15) .P 9 (5) .P 9 (2) .P( 2) , 

9 8 = 28150- Pj M ^P«.P< a \ 

97 = 2111 • (at + 16608) • P 7 (64) • P 7 (2) , 

9 6 = 21300 • + 15054) • (a: + 20971) • P 6 (37) • P 6 (23) • P 6 (2) • pf } , 

95 = 8699 • (x + 1134) • P 5 (57) • P 5 (4) • P 5 (3) , 

94 = 7621 • (x + 27997) • P 4 (7) • P 4 (42) • P 4 (8) • P 4 (6) , 

93 = 18283 • + 26460) • P 3 (58) • P 3 (4) , 

92 = 2235 • (at + 8688) • (x + 20285) 

x (x + 24023) • P 2 (34) ■ P 2 (12) ■ P 2 (8) • P 2 (3) ■ P 2 (2) , 
9i = 2139 • {x + 19284) • (a: + 19339) • P^ 55) ■ P^\ 
9o = 23255 • (x + 30075) ■ {x + 13139) • P (41) ■ P (15) ■ P (2) , 

where the polynomials Pm^ or iftf are polynomials of degree N : 

PW = ^ + ... (A.4) 

Appendix B. Singular behavior of Xd^ an d Xd 

Using the matrix connection method [5] , we have obtained the values of the amplitudes 
at the singularities, lying on the unit circle |i| = 1, of pQ and X^P ■ 

(3) 

The linear differential equation for \ d has the singularities x = 0, 1, — 1, do and 
the roots of 1 + x + x 2 (t — x 2 ). The singular behavior of x^ a t the singular point 
x s , denoted by x\ (singular, x s ) read (in the local variable u = x — x s ): 

Xf (singular,!) = Q + J- + ■ \ + ~ M«) (B.l) 
with a = 0.469629259 • • • 



*? (singular,-!) = ^ ln(„) 2 + - ii^l ] - ln( « ,. 



X^ 3) (singular, x ) = -i (1 - j) (1 + i y/l) ■ b ■ u 7/2 
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where x = —1/2 + z y/3/2 and with b = 0.203122784 • • • 

The linear differential equation for xP has the singularities t = 0,1,— l,oo. 
Denoting by it = t — t s the local expansion variable, where t s is the singularity, 

Xf (singular, 1) = (-L - i - |V i 



the singular behavior for x'P read 



8tt 



1 . , . 2 8 ln(2) - 7 -i2n , . . 
ln(u) 2 + V ; o ln(it) 



with c = -1.120950429 • • • and: 



xP (singular,-!) = —^—^. u 7 -\n{u). 
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